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Abstract 

We compute the full one-loop EW contributions of 0{a^Q^y^ entering the electron- 
positron into two 6-quarks and one gluon cross section at the Z peak and LC ener- 
gies. We include both factorisable and non-factorisable virtual corrections, photon 
bremsstrahlung but not the real emission ofiy^andZ bosons. Their importance for 
the measurement of as from jet rates and shape variables is explained qualitatively 
and illustrated quantitatively. Their impact on the forward-backward asymmetry 
is also analysed. 

1 Introduction 



Jet samples emiched in 6-quarks produced in e^e~ annihilations are used for sophisticated 
tests of QCD, primarily because they enable one to distinguish between quarks and gluons, 
thanks to 6-flavour tagging (typically, by exploiting high px leptons and/or microvertex 
techniques). This is unlike the case of lighter flavour^. To name but a few examples, by 
studying 6-jet samples, one can: (i) verify the flavour independence of as; (ii) study the 

^Work supported in part by the U.K. Science and Technology Facilities Council (STFC), by the 
European Union (EU) under contract MRTN-CT-2006-035505 (HEPTOOLS FP6 RTN) and by the ItaUan 
Ministero dell'Istruzione, deU'Universita e della Ricerca (MIUR) under contract 2006020509_004. 

^With the possible exception of c-quarks, whose tagging efficiency is however much lower in comparison 
to that of 6-quarks, so as to make the former much less suitable than the latter to phenomenological 
investigation. 
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properties of the QCD force carrier (the gluon); (iii) measure the 6-quark mass. (For an 
authoritative review on experimental tests of QCD in e~^e~ events see [I] and references 
therein.) 

Of particular relevance is the three-jet sample, as it impinges on all types of analysis 
(i)-(iii). It is therefore of paramount importance to give predictions for the e~^e~ — > 
bbg cross section to the highest degree of precision, thereby necessarily involving the 
computation of all higher order corrections within the Standard Model (SM). Whilst the 
computation of QCD effects through one-loop has been tackled some time ago [21 E] the 
case of Electro- Weak (EW) corrections is not available in the literature. We remedy this 
shortcoming here, by calculating the full one-loop EW corrections to bbg observables in 
electron-positron annihilations generated via the interference of the graphs in Figs. 1-6 
of Ref. [1] (see also [5l [6]) with the tree-level ones for e~^e~ — 7*, Z — bbg. In doing so, 
notice that we will be including photon bremsstrahlung; in contrast, we will refrain from 
computing real and Z boson radiation, as we will argue that this may not enter the 
experimental jet samples. 

Finally, notice that, while QCD corrections are dominant at low e'^e~ energy, EW 
ones become relatively more and more important as the latter grows larger, because of 
surviving Sudakov logarithms from which QCD interactions are immune. Besides, EW 
corrections also carry the hallmark of parity- violating effects, which are generally peculiar 
to new physics beyond the SM, so that they ought to be accounted for in its quest. 

The plan of the rest of the paper is as follows. In the next Section, we describe the 
calculation and give the input parameters and jet-selection algorithms. Then, in Sect. [HI 
we present our numerical results for LEPl/SLC, LEP2 and a future Linear Collider (LC). 
We conclude in Sect. HI 

2 Calculation 

The procedures adopted to carry out our computation have been described in to 
which we refer the reader for the most technical details. Here, we would only like to point 
out that we have neglected the masses of the fe-quarks throughout. However, whenever 
there is a boson in the virtual loops, account has to be taken of the mass of the top 
(anti)quark, which we have done. Furthermore, before proceeding to show our results, 
we should mention the numerical parameters used for our simulations. We have taken 
the top (anti)quark to have a mass rrit = 171.6 GeV. The Z mass used was Mz = 91.18 
GeV and was related to the mass, Mw, via the SM formula = Mz cos 9w, where 
sin^ 6\Y = 0.222478. The Z width was Tz = 2.5 GeV. Also notice that, where relevant, 
Higgs contributions were included with Mh = 115 GeV. For the strong coupling constant, 
as, we have used the two-loop expression with Aq,^^^'' = 0.325 GeV in the MS scheme, 
yielding af^{Ml) = 0.118. 
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As for the jet definition, partonic momenta are clustered into jets according to the 
Cambridge jet algorithm [7] (e.g., when i/ij < i/cut with ycut = 0.001), the jets are required 
to he in the central detector region 30° < ^jcts < 150° and we request that the invariant 
mass of the jet system Mf^g is larger than 0.75 x -^i. If a real photon is present in the 
final state, it is clustered according to the same algorithm, but we require that at least 
three "hadronic" jets are left at the end (i.e., events in which the photon is resolved are 
rejectedjfl. We further make the assumption that both 6-jets can be tagged (as described), 
including their charge (e.g., via the emerging lepton or the jet-charge method). (For sake 
of illustration we take the efficiency to be one.) In order to show the behaviour of the 
EW corrections we are calculating, other than scanning in the collider energy, we have 
considered here the three discrete values of ^/s = Mz (also in view of a GigaZ option of 
a future LC), ^/s = 350 (the tt threshold) GeV and y/s = 1 TeV (as representative of the 
Sudakov regime). 



3 Numerical Results 

In Fig. [H we present the effects on the cross section (integrated within the experimental 
cuts defined in the previous Section) induced by different terms of the order asa^u con- 
tribution relative to the lowest-order cross section for e+e~ bbg, plotted as a function 
of the Centre-of-Mass (CM) energy, in the range from 150 GeV to 1 TeV. The curves 
represent the effects of the QED (virtual and real) corrections only, the gauge bosons 
self-energy corrections, the non-factorisable graphs involving four- and five-point func- 
tions with WW exchang^, the weak corrections with the non-factorizing WW graphs 
removed and the sum of the previous ones. Notice that the total effect is increasingly 
negative, as ^/s gets larger, reaching the —13% or so level at 1 TeV. As already stressed 
in Ref. [Ij and visible from the plot here, such a big negative correction is mainly due to 
the presence of the WW non-factorisable graphs, which develop the aforementioned large 
Sudakov double logarithms in the very high energy regime. The pattern of the various 
corrections seen here is not very different from the one seen in Ref. [1| (Fig. 7 therein) 
with the notable difference that in the case of 6-samples one can appreciate the onset 
of the ti virtual threshold, which was instead invisible in the fully flavoured sample of 
Ref. m. 

The ability to efficiently tag 6-quark jets enables one to deflne observables in bbg flnal 
states which are not (easily) reconstructable in the case of the full three-jet sample. One 

•^As explained in [3], this serves a twofold purpose. On the one hand, from the experimental viewpoint, 
a resolved (energetic and isolated) single photon is never treated as a jet. On the other hand, from a 
theoretical viewpoint, this enables us to remove divergent contributions appearing whenever an unresolved 
gluon is ejected via an infrared (soft and/or coUinear) emission, as we are not computing here 0{asa'^yi) 
one-loop QCD contributions to e+e~ — > bbj. 

''This is a gauge invariant subset of the complete corrections. 
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Figure 1: Relative effect on the integrated cross section due to different contributions to 
the order a = oem correction, as a function of the CM energy. 
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Figure 2: bb invariant mass distribution at the Z peak. 
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Figure 3: bb invariant mass distribution at 350 GeV. 




Figure 4: bb invariant mass distribution at 1 TeV. 
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Figure 5: (1 — T)^ distribution at the Z peak. 



example is the invariant mass of the bb pair, Mf,i, which we plot in Figs. [2H11 Here, the 
largest contribution to the total correction comes from QED Initial State Radiation (ISR), 
primarily because of the radiative return phenomenon. 

In fact, the correction is up to 20% for very high energies and intermediate invariant 
mass and -30% for large invariant mass whereas it is a (constant) —45% at small energy. 
The purely weak corrections are negative and at the level of —15% at most, for all en- 
ergies. The higher order QED radiation tends to compensate the order oem effect, since 
it enhances the cross section. One should also notice that, again thanks to a Sudakov 
effect, by raising the CM energy, the relative weight of the weak corrections becomes 
more important and the effect of higher order QED corrections diminishes. Finally, it is 
worth noticing that, far from the Z peak, the cut M^j^ > 0.75 x ^/s is more effective in 
suppressing the radiative return phenomenon, reducing in turn the relative effect of ISR. 

In the following some event shape variables are considerec^]: the thrust T [11] and the 
C-parameter [12] (see Ref. [13] for their definitions). In Figs. [5l [6] and [TJ the spectrum 
of (1 — T)^ is shown. This distribution is one of the key observables used for the 
measurement of as in e~^e~ collisions [13]. It is worth noticing that whilst the purely 
weak corrections give an almost constant effect on the whole T range, the presence of the 

^Recalling footnote 3, one should notice that these observables will depend on ycut- However, this 
reflects standard experimental procedures [10] aiming at removing resolved photons from the hadronic 
sample and we also have verified that the relative size of the EW effects computed here does not depend 
on ycut- 
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Figure 6: (1 - T)f^ distribution at 350 GeV. 
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Figure 7: (1 - T)^ distribution at 1 TeV. 
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Figure 8: ^ distribution at the Z peak. 



real bremsstrahlung gives a non-trivial effect in the region T > 0.95 (at least at higher 
energies). The numerical impact of the various contributions is similar in the case of the 
C-parameter. Here, the region where QED real radiation introduces non-trivial effects is 
C < 0.1. Given the size of the various EW corrections (the QED ones being of tens of 
percent), including the purely weak ones (which are steadily at -7% or so) and in view of 
a precise measurement of as from bbg samples at future LCs, it is clear that such effects 
will play an important role and thus cannot be neglected. 

In Figs. [11], [12] and [13] the Cambridge y distribution is shown. For each event, the 
observable y is defined as the minimum (Cambridge) yij such as yij > ycut = 0.001. 
Also on this distribution the weak effects are quite constant and between -5% and -10% 
(increasing with energy), while real radiation effects are significantly larger (and negative) 
and distort the LO shape. 

Figs. [131 [15] and [16] present instead the cross sections integrated over y in the range 
ycut < y < ymax, as a function of ymax- Corrections can be very large in such distributions, 
generally at any energy, reaching the minus several tens of percent (QED ones) or nearly 
the -10% (weak ones) level. 

If one combines fe-(anti)quark flavour tagging with jet-charge measurements, it is pos- 
sible to define the forward-backward asymmetry, ApB, for the case of 6-jet samples. As a 
clear anomaly in this observable (defined already at LO in the case of bb final states) has 
survived after the LEP and SLC era [TU [TJ], it is worthwhile to investigate the impact 
that (hitherto neglected) EW contributions through 0(asa|]yi) can potentially have in 
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Figure 9: % distribution at 350 GeV. 
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Figure 10: ^ distribution at 1 TeV. 
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addressing the discrepancy between data and SM. Here, we use the following definition 
for the asymmetry distribution as a function of the bb invariant mass: 

r^/"^ jn da rTT jn da 

AFB[Mf,i) - da 

JO "'^^de^dM^i 

with 6}, being the polar angle of the 6-jet relative to the electron, wherein both the numer- 
ator and denominator (normalisation) involve only terms of 0{aa,oi\y^ plus Oi^a^a^y^) 
and where the final state comprises three jets. The integrals involved in Eq. ([1]) are over 
the phase space allowed by the cuts defined in the previous Section. Whilst the size of 
the asymmetry is substantial through the computed orders, and the full 0(as«EM) rates 
are sizeably different from the 0(q;sQ^em) ones (at least at very small and very large ener- 
gies), our results should clearly be folded with all known higher order corrections to the 
e^e~ —>■ bb process (work on this is in progress). 

Finally, in Fig. [201 the forward-backward asymmetry, now integrated over M;,^, 

r'"^/^ J/i da rn in da 
Jo ^^b^^ 

is plotted as a function of the CM energy and as obtained at lowest order and by including 
only QED corrections or the full one-loop EW corrections. As the CM energy raises, the 
increasing effect of the purely weak corrections is more and more evident. 



4 Conclusions 

A careful analysis of actual e~^e~ —>■ bbg data is in order then, involving one- loop EW 
effects. In this regard though, one caveat should be borne in mind: as emphasised in 
Refs. [H [9] (albeit in the hadronic context), particular care should be devoted to the 
treatment of real and Z radiation and decay in the definition of the jet sample, as 
this will determine whether tree-level and Z bremsstrahlung effects (neglected here) 
have to be included in the theoretical predictions through 0(as«EM)- However, given 
the cleanliness of jet samples produced in electron-positron machines, as compared to 
hadronic ones, we believe that the former contribution can effectively be disentangled, 
thereby rendering our present predictions of immediate experimental relevance. 
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Figure 11: Cambridge y distribution at the Z peak. 
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Figure 12: Cambridge y distribution at 350 GeV. 
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Figure 13: Cambridge y distribution at 1 TeV. 




Figure 14: Cross section as a function of the Cambridge maximum y at the Z peak. 
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Figure 15: Cross section as a function of the Cambridge maximum y at 350 GeV. 




Figure 16: Cross section as a function of the Cambridge maximum y at 1 TeV. 
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Figure 17: Aps distribution at the Z peak. 




Figure 18: ApB distribution at 350 GeV. 
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Figure 19: ApB distribution at 1 TeV. 
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Figure 20: ApB as a function of the CM energy, obtained at tree level, adding only QED 
corrections and the full one-loop EW corrections. 
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